Spinning compact binary inspiral II: Conservative angular dynamics 

Laszlo Arpad Gergely^'^* 

^Department of Theoretical Physics, University of Szeged, 

Tisza Lajos krt 84-86, Szeged 6720, Hungary 

^Department of Experimental Physics, 

University of Szeged, Dom ter 9, Szeged 6720, Hungary 

* E-mail: gergely@physx.u-szeged.hu 

(Dated: November 23, 2010) 

We establish the evolution equations of the set of independent variables characterizing the 2PN 

rigorous conservative dynamics of a spinning compact binary, with the inclusion of the leading 

order spin-orbit, spin-spin and mass quadrupole - mass monopole effects, for generic (noncircular, 

nonspherical) orbits. More specifically, we give a closed system of first order ordinary differential 

equations for the orbital elements of the osculating ellipse and for the angles characterizing the spin 

orientations with respect to the osculating orbit. 

("^ ' We also prove that (i) the relative angle of the spins stays constant for equal mass black holes, 

Cn ' irrespective of their orientation, and (ii) the special configuration of equal mass black holes with 

►^ , equal, but antialigned spins, both laying in the plane of motion (leading to the largest recoil found in 

Q ' numerical simulations) is preserved at 2PN level of accuracy, with leading order spin-orbit, spin-spin 

h^ ' and mass quadrupolar contributions included. 

o 

CN ; I. INTRODUCTION 

O ', Compact binaries composed of neutron stars or stellar size black holes are among the most likely sources to emit 

^ ■ gravitational waves in the frequency range of the Earth-based gravitational wave detectors LIGO and Virgo |l|. 
$_i ' Supermassive black holes in the mass range of 3 x 10^ ^ 3 x 10^ solar masses reside in the centers of galaxies and 
!^i: following the merger of their host galaxies, they also merge. In the process they create powerful gravitational waves, 

detectable in the lower mass range by the space mission LISA 2] . 

pO ' By definition the inspiral is the regime of the orbital evolution, during which the post-Newtonian (PN) parameter 

^ ' e = Gm/c^r ss (v/c) (where m = mi -\- 7TI2 is the total mass, r and v the orbital separation and relative velocity 

^^ of the binary) is small and where the leading order dissipation is due to the gravitational waves. As two galaxies 

^T^ ' merge, their supermassive black holes are subject to both gravitational radiation and dynamical friction. The former 

overcomes the latter at about ei„ = 10~^ (the actual number only weakly depends on both the stellar distribution and 

mass Q). During the inspiral which follows, the parameter e increases. When e approaches its value at the innermost 

stable orbit, the PN description becomes increasingly nonaccurate, therefore the subsequent plunge is better described 

by numerical evolutions, or as an alternative, by expressions traced back to the PN approach, arising either from the 

effective one-body model, calibrated to numerical relativity simulations [J| or from a phenomenological transition 

phase, with coefficients again calibrated by comparison with specific, numerically generated waveforms [5|. Finally, 

the ringdown follows, when the newly formed black hole radiates away its physical characteristics, with the exception 

of mass, spin and possibly electric charge (for a review of quasinormal modes of black holes see Ref. [y]). 

\^ • The spin and quadrupole moment of the supermassive black hole at the Galactic center can be measured via 

5^ \ astrometric monitoring of stars orbiting at milliparsec distances [7] , and this can also be a test of the general relativistic 

no-hair theorem. 

The spin affects the horizon of the black holes, therefore those observations on black holes which indicate the size of 
the horizon, will also lead to indirect spin magnitude estimates. (Estimating the quadrupole moment from observing 
a two-dimensional projection of the horizon would be less straigtforward.) Both stellar size and supermassive black 
holes can have accretion disks and jets in their active periods. Whenever these observations are connected to the 
presence of a jet and the direction of it can be identified (projection effects may again obstruct this), we also obtain 
information on the black hole spin axis, assuming it is aligned with with the symmetry axis of the magnetic field and 
hence the jet direction. Jets from rotating black holes have been shown to be stable jSjj. Spin direction can be also 
inferred from observations on the radiation of the accretion disk. 

Such observational spin estimates can be made at least by four methods: 

i) Reverberation mapping of the observed optical/X-ray lines (highly excited Mg, O, C) in active galaxies to 
determine the radius and velocity pattern of the Broad Line Region. This depends on the metric, decreasing with 
increasing spin. From such considerations the mass, spin and spin orientation of the black holes can be estimated [9|. 
In particular, information on the spin direction of the central black hole of the Seyfert galaxy Mrk 110 was obtained 
by estimating the central black hole mass in two independent ways. First, assuming that the broad emission lines are 
generated in gas clouds orbiting within an accretion disk, the mass could be determined as function of the inclination 
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angle of the accretion disk. Secondly, detecting the gravitational redshifted emission in the variable fraction of the 
broad emission lines, a central black hole mass, which is independent on the orientation of the accretion disk could 
be deduced [U^ . 

ii) There is a strong effort towards imaging with millimeter Very Long Baseline Interferometry (VLBI) the event 
horizons of Sagittarius A* (SgrA*) and Virgo A (M87), which again depend on the spin. For SgrA*, the radio source 
at the center of our Milky Way, millimeter and infrared observations require the existence of a horizon [ll| . Analysing 
the peaks of the power density spectra in the light curves of X-ray flares from the Galactic Center black hole the 
mass and spin were inferred [12i]. A compact emission region (bright spot) in a circular orbit and the lightcurves of 
its observed flux and polarisation depend on the mass and s pin of the black hole. The emitted polarisation fraction 
is polarised orthogonally to the spin axis of the black hole [Ij]. Unlike SgrA*, M87 exhibits a powerful radio jet, 
allowing future VLBI data to constrain the size of the jet footprint, the jet coUimation rate, and the black hole spin 

iii) The size of the jet launching region in Active Galactic Nuclei (AGN) is determined by the Blandford-Znajek 
effect, which in turn depends on the spin [l^- Measurement of the diameter of the jet base (e.g. in M87) gives 
evidence for small sizes, regarded as signature of a large spin |16l |. 

iv) The low-energy cutoff in the energetic electron spectra of the jets suggested by the radio spectra [l3| is con- 
veniently explained by the pion decay resulting from proton-proton collisions [181 . The latter mechanism require 
a relativistic temperature in the accretion disk near the foot of the jet, which translates to the central black hole 
spinning extremely fast jl9| . 

Some of these methods will certainly also work for stellar size black holes due to scale invariance arguments in 
accretion phenomena. The jet/disk geometry has been constrained for the stellar black holes XTE J1118+480 and 
GX 339-4 120]. 

From all these observations we conclude that it is necessary to include the spin and quadrupole moment of black 
holes when modeling their binary systems. 

In this paper we investigate the 2PN rigorous conservative dynamics during the inspiral of a spinning compact binary 
system, by including leading order spin-orbit (SO), spin-spin (SS) and mass quadrupole - mass monopole (QM) effects, 
for generic (noncircular, nonspherical) orbits. Because of these interactions the spins undergo a precessional motion 
[2l|-[23- Various aspects related to the leading order contribution to both the conservative and dissipative part of 
the dynamics due to the SO interaction were discussed in Refs. I23| -f25i. while the corrections represented by the 
SS coupling in Refs. [23-[2al, and by the QM coupling in Refs. [23|-[29]. The radial motion under the Newton- 
Wigner-Pryce spin supplementary condition [30|, with all these contributions included is fully solved in Ref. [3l|. 
The Hamiltonian approach including spins has been also widely discussed |32| |. Based on numerical work, empirical 
formulas for the final spin have been proposed in Refs. [33|. Zoom- whirl orbits, which were known to exist for particles 
orbiting Kerr black holes [3J] , also appear in the framework of the PN formalism [3a| , their likeliness increasing with 
the spin [3^]. Gravitational wave emission can lead to a spin- flip [33]- [331 in X-shaped radio galaxies [33], [S^j- It is 
a more recent result that for mass ratios 0.3 -^ 0.03 the combined effect of SO precession and gravitational radiation 
will result in a substantial spin-flip already in the inspiral phase [3|, [40|. 

With the spins and mass quadrupole moments included, the number of variables in the configuration space increases 
drastically. In Ref. [4l| a minimal and conveniently chosen set of independent variables for such a system was 



established. Notational correspondence of some of these variables with quantities employed in Refs. [42|-[4J], where 
the dynamics of spinning compact binaries (without mass quadrupolar contributions) has been also discussed, is 
established in Appendix [Xj Beside the masses mi, the number of independent variables characterizing the total and 
orbital angular momenta (J and L, respectively) and spins Si was shown to be 6, chosen either as 

(a) 3 angles (spanned by the Newtonian orbital angular momentum Ln with the total angular momentum J 
and with the spins Si, denoted as a and k^, respectively) and 3 scales (the normalized magnitudes of the spins 
Xi = (c/G) (Si/mf) and the magnitude of the total angular momentum J), or equivalently as 

(b) 5 angles and a scale. In this case the dimensionless spin magnitudes Xi could be replaced by the azimuthal 
angles ipi of the spins, measured in the plane of motion from a suitably defined node line 1 (the intersection of the 
planes perpendicular to the total orbital momentum J = L-I-S1-I-S2 and to the Newtonian orbital angular momentum 



Ln). The relation between the two sets of variables is given by Eqs. (46)-(47) of [41]. 

In the present paper we discuss the conservative dynamics of these two sets of independent variables. All lengths 
involved (J and Xi) are constants of motion, as J is conserved to 2PN accuracy [23] and the spins undergo a precessional 
motion [2l|-(22lj. Therefore our goal reduces to the study of the dynamics of the angular variables. In the process we 
also derive the evolutions of the parameters (a^, e^) of the osculating ellipse; of the spin relative angle 7 span by Si 
and S2; and of the periastron, given by the Laplace- Runge-Lenz vector An. 

We start with a discussion of the evolutions under a generic perturbing force in Sec. [TTl First we monitor how the 
Keplerian dynamical constants evolve. This allows us to determine both the evolution of (a^, e^) and of the vectors 
Ln, A-n, and ascending node vector 1. We also determine here how the evolution of the true anomaly parameter 



Xp (measured from An to the actual location f of the reduced mass particle) is modified by the perturbing force. 
The specific perturbing force components generated by PN, 2PN, SO, SS and QM effects are listed in Appendix iBl 
together with the components of the spin precession angular velocity. 

Employing these results, also the spin evolution equations discussed in detail in [41j, we are able to derive in Sec. 
mil the evolution of a. Equation (14) and (15) of 41] show that once the evolution of Xp and a are established, the 
evolution of the angle ipp measured from 1 to An, and of the angle —4>n measured from 1 to an arbitrary inertial axis 
X ± J (see Fig 1 of [41]) also follow, which complete the characterization of the evolution of the Euler angles. Then, 
in Sec. IIVI we derive the evolutions of Ki, 7, and ipi. With this we fulfill the task of characterizing the evolution of 
the variables composing the independent sets (a) and (b). 

We discuss special spin configurations in Sec. FVland present our Concluding Remarks in Sec. I VII 

Notations and conventions. The gravitational constant G and speed of light c are kept in all expressions. For any 
vector V we denote its magnitude by V and its direction by V. 

The reduced mass is /i = mini2/m. We assume that mi > TO2, thus the mass ratio v = 7712/7711 < 1 and the 
symmetric mass ratio rj = fi/m — v / {\ A- v) G [0, 0.25]. 

The mass quadrupole moment originates entirely from rotation, being therefore characterized by a single quadrupole- 
moment scalar Qi = — {CP /c^^ wxirn^, with the parameter w € (4, 8) for neutron stars, depending on their equation 
of state, stiffer equations of state giving larger values of w [23], [43| and w = 1 for rotating black holes ^. The 
negative sign arises because the rotating compact object is centrifugally flattened, becoming an oblate spheroid. 

The inertial system /Ci has the arbitrary inertial x axis x and J as its z axis. We also define the noninertial systems 
/Cl and ICa with Ln as the common z axis, the x axes being 1 and An, respectively. Then the y axes are m = Ln x 1 
for /Cl and Qn = Ln x An for ICa- 

II. EVOLUTIONS IN TERMS OF A GENERIC PERTURBING FORCE 

Although there is no notion of gravitational force within general relativity, in the PN regime the motion of a 
compact binary can be regarded as a perturbed Keplerian motion, with perturbations coming from the difference in 
the predictions of general relativity with respect to Newtonian gravity. Therefore one can adopt the terminology of 
celestial mechanics, regarding the modifications induced by general relativity as perturbing forces. 

Any perturbed Keplerian motion is characterized by an acceleration 

Gm^ , , . 

a= ^r + Aa. (1) 

We find convenient to express Aa in the basis ICa with basis vectors {f(i)} — (An, Qn, Ln) as 

3 
Aa=^a,f(i). (2) 



A. Keplerian dynamical constants 

Starting from the definitions of the Keplerian constants of motion E^^ = fiv'^ /2 — Gm^/r, Ln = /ir x v, and 
An = V X Ln — Gmiir, it is straightforward to show that 

i;jv = Aiv • Aa , (3) 

Ln = MI" X Aa , (4) 

An — Aa x Ln + v x Ln 

= /i [2 (v-Aa) r- (r-Aa) V- (r • v) Aa] . (5) 



By employing the decomposition of r and v in the basis ICa, given by Eqs. (|Bl[) - (|B2p . also the decomposition of the 
perturbing force acting on the unit mass ^, and finally the generic formula for the time-derivative of any vector V, 

V ^VV+V^V , (6) 



we obtain for the magnitudes 



and for the directions 
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Gm^^ ^ {A]sj + Gmfj, cos x„) 
-ai — smxp + a2 



Ln = (a2Cosxp -aisinxp)^^ , 

• r , f ■ . M''(^JV + Gm^cosXp) 

An = a2LN + [a2C0SXp- aismxp) 

Ln 



a3-—r (sin XpA.Tss -cos XpClTss) , 
Ln ^ ' 

Ln , Gmy.^ . . ' 

-fli -; h ^— r sm Xp (02 cos Xp - ai sm Xp) 



QN-as-j—smXpLTss , 
Ln 



' An LnAn 
where r is given in terms of the true anomaly parameter Xp by the standard formula (jB20[) 



(7) 



(8) 



B. Radial semimajor axis a^ and radial eccentricity e, 

We note that the constraint A\, = (Gnifi) +2EL%/fi is preserved by the evolutions ([7]), therefore only two of these 
equations are independent. From them we can also derive evolution equations for the parameter pr = Lj^/Gm^^ 
and eccentricity Cr = An /Gm^ of the conic orbit. For bounded orbits we could introduce the semimajor axis a^ — 
Pr/ (1 — e^) = L'jf/GmiJ^ (l — e^) = —Gmfi/2EN of the osculating ellipse instead, and derive evolution equations for 
the pair (a^, e^). In this way we obtain two Lagrange planetary equations: 
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Or (1 - e^) 



Gm 



-aisinxp + a2 (er + cosXp)] , 



02 (1 + 2er- cos Xp + cos^ Xp) - »! (er + cos Xp) sin Xp 
(1 + erCOSXp) 



(9) 



(10) 



Here we have employed the true anomaly parametrization (|B20p written in terms of osculating ellipse orbital elements 



Qr (1 - e^) 

1 + Cr COS Xp 



(11) 



C. The noninertial system Ka 

Rewriting Eqs. (|8|) in the form of precession equations by inserting An = Qn x Ln, Qn — —An x Ln in the first 
expression and Qn = Ln x An, Ln — — Qn x An in the second; also computing the time derivative of Qn from its 
definition gives 



f(i) = ^A X f(i 



(i) 



with the angular velocity vector 
[ir cos Xp 



n. 



a-A- 



^N 



- fir sm Xp /^ 
An + as — Qn 



^N 



ai 



A 



N 



(aismXp-a2C0SXp, 



Gmp?r sin Xp 
LnAn 



-IN 



(12) 



(13) 



With this we have established the time evolution of the noninertial basis K-a- 

The PN order of Jl^ is O {^a) — e^^l'^O {ai/cj. Employing the contributions to a^ from Appendix [B] and Eq. (58) 



of |4li one finds 






o (f^r^) 



O(niO) - o(e3/2 



E^(l'- 



2/£-3^ 



Xk 



Lfe=l 



o(r~i) 



0(f2f) = 0{e')0{r^)xiX20{T-^) 



{nf^) ^ O {e') O (r;) 



cm 



■ 2 

E' 

.fe=i 



(.^'^-3) ^,x^ 



©(T-i) 



(14) 



with r being the radial period, defined as twice the time elapsed between consecutive f — Q configurations. 
A couple of immediate remarks are in order: 

(1) If 03 = (no perturbing force is pointing outside the plane of motion), Ln (the plane of motion) is conserved, 
while both An and Qn undergo a precessional motion about Ln (in the conserved plane of motion). 

(2) If oi = a2 = (the perturbing force is perpendicular to the plane of motion), then An undergoes a precessional 
motion about Qn and vice-versa, while Ln precesses about r. 

D. True anomaly Xv 

As the basis{f(i)} is comoving with the plane of motion and the periastron, the position vector r — a;*f(i) [with x^ 

given by Eq. (|B1[) ] changes according to v = i;*f(i) + a;*f(i) = 2;'f(i) + x^VIa x f(i). A straightforward computation, 
employing 



then leads to 



From here 



iXp-J'XpSmXp , 



^N = /"?■ 



Xp 



X = f sin Xp + rxp cos Xp 
(flA ■ Ln 



0, 



^N 



Xp 



IftA ■ LnJ 



Ln 



(15) 
(16) 

(17) 



Therefore the deviation from the Newtonian expression is due to the component of il^ along Ln- The importance 
of Eq. (|17p lies in allowing to pass from time derivatives to derivatives with respect to Xp in the evolution equations 
(HI), ((9|)- (fT0)) . which then become ordinary differential equations. 
It is also immediate to derive w^ and calculate E^ as 



En = 7^ h fJ-r Xp V^A ■ Ln 



2 r 

By inserting Eq. (J17p . we obtain the radial equation 



fir 



IUa ■ Ln j 



2En 2Gm L% 



^J■ 



9 9 



(18) 



(19) 



Remarkably, all terms arising from the precession of the basis vectors cancelled out and we formally recovered the 
radial equation for the Keplerian motion. This is not surprising, as the dynamical quantities En, Ln refer to the 
osculating ellipse. 



E. Ascending node 1 

The basis vectors of JCl are related to the basis vectors of ICa by a rotation in the x-y plane with angle — V'p, thus 

1 = cost/i^An — siuT/ipQN , (20) 

rh = sinV'pAN + cost/ipQn • (21) 



The time derivative of the direction of the ascending node is therefore found as 

— 1 = -V'p (sinV'pAN + cos'(/'pQnJ +cosV'p— An - sinT/ip — Qn 
= cos ipp ( flA — ^pLn ) X An — sin ipp ( flA — tppi^N ) x Qn 

= (Qa - i^pLn) X i . (22) 

Similarly we can derive the evolution of m as 

— m = ^pp f cos V'pAn - sinV'pQNJ + sinf/'p— An + cos?/'p-t-Qn 
= cos Ipp ( flA — i/'pf'N j X Qn + sin -^p ( ft a — V'pl'N ) x An 
= f ri^ - ''/'pLn) X m . (23) 

As it was to be expected, the unit vectors 1 and rh undergo a precession characterized by the angular velocity vector 

nL = nA~ v^pLn . (24) 



III. EULER ANGLE EVOLUTIONS 



Now we have all necessary elements for deriving the evolution of the angles which enter the set of indepen- 
dent variables. First we remark, that the time derivative of the definition of the argument of the periastron 

^p — arccos (1 • An ) , by employing Eqs. (fT2)) and p2|) gives an identity. 



A. Inclination a 



From the definition of the inclination a = arccos I J • Ln ) , employing the constancy of J up to 2PN [23| and the 
derived precession equation for Ln we find 

— sina a — J-— Ln = J- ( ^A x Ln ) = ^A' (Ln x J J = — sina Sl^-l , (25) 

thus 

^^../rcosj^bp + xr,) ^ (26) 

Ln 

B. Longitude of the ascending node —(j>„ 

By employing Eq. (14) of ^l[ we find the evolution of the azimuthal angle — </>„ of the ascending node 1 as 

• iirsinj^p + Xp) .„„^ 

(Pn = -0.3 r ■ ■ (27) 

Ljv sm a 

Quite naturally, both the orbital inclination and the ascending node can be changed only by a force perpendicular to 
the orbit. 



C. Argument of the periastron t/j. 

From Eq. (15) of [41| and Eq. (I?fl) the evolution of ipp + Xp emerges as 



• Ln fir sin {ipp + Xp) .^„. 

^P + Xp = ^-a3 ^—r —■ (28) 

fir'^ Ln tan a 



Again, only the perturbing force component along Ln contributes. Combining Eqs. ((28)) and (|17p leads to the 
evolution equation of the third Euler angle. 



J^A • Ln - ^p = 03 



fir sin {-tpp + Xp) 
Ljv tan a 



(29) 



The left-hand side is 17^, • Ln, such that the unit vectors 1 and ih undergo a precession characterized by the angular 
velocity vector 



Hl = as 



fir 



A , ■ A , sin(V'p + Xp)f 
cos XpAn + sm XpQn H -^^ — ^^n 



tana 



(30) 



The first two terms of the bracket combine to f . If there is no perturbing force perpendicular to the orbit, 1 and m 
stay unchanged. 

The evolution of i/'p in detail reads 



• Ln , . Gmfi^rsinxp ^r sin (■0p + Xp) 

iPp = -fli-i (aismXp-a2COSXp) z — -. 03 — 

An LnAn Ln ta-na 



(31) 



Equations (^5]) . (I?fl) and (PT|) are Lagrange planetary equations for the angular orbital elements. With the use 
of Eq. (|17p . by passing from time derivatives to derivatives with respect to Xp: these become ordinary differential 
equations. During the inspiral the perturbing force components a^ arise as a combination of relativistic (PN and 
2PN), SO, SS and QM contributions, and are given in Appendix [Bl 



IV. SPIN ANGLE EVOLUTIONS 



A. Spin polar angles Ki 



The spin polar angles Ki = arccos ( Si • Ln ) evolve due to the spin precessions (see Appendix |B]) and the evolution 
of Ln, as 

(32) 



(33) 



Q,A X Ln j -Si + Ln- [p.-, x S;) = [VIa - rj;) • ( Ln x S; 
In order to proceed, we need the expression (|B3p of the spin, such that 



Ln x Si = sinKi 



sin (V'p - V'i) An + cos (V'p - V'O Q 



N 



and we find 



ki = [Vli ■ An j sin (^p - ipi) + f Jli • Qn j cos (^p - ipi) - a^— sin (^p + Xp - i^i) 



(34) 



The relative orientation of spins with respect to the orbital angular momentum is unchanged only if the perturbing 
force lies in the plane of motion (03 — 0) and if the spin precession axis is along Ln- The latter condition is obeyed 
by the SO precession, but not by its SS and QM corrections (except for perfect perpendicularity of the spins to the 
orbital plane, when also 03 = holds, see Appendix iBl thus ki — 0). Starting from this and the remark 03 oc C (e^^^), 
also the estimates (JB35I) we find 



O 



(k.) = O (e^/^) O (rf) hx. + O (z.2-3) ^^ ] ^(^-1) . 



(35) 



B. Relative spin angle 7 



For this we take the derivative of its definition 7 = arccos (Si ■ S2 ) and obtain 



sin 7 7 — (Oi — O2) • (Si X S2 



(36) 
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If the spins are either ahgned or antiahgned with each other, such that Si x S2 = 0, then 7 = 0, irrespective of the 
mass ratio. 

Otherwise, by employing Eqs. (56) of 41] and also (Si x S2) • Sj = 0, we rewrite the condition ([55t as 



3G 



sin7 7 - ( ^^ ^ ^ Ln + {f ■ [{l-W2V-^) S2 - (1 - wiv) Si] } f ) • (Si x 83 



Equal mass {v — 1) black holes {wi = 1) trivially imply 7 = 0, irrespective of the orientations of the spins. 
For the generic case from Eq. (jB3p we have 



Si X S2 = [cos Ki sin K2 sin ('0p — "02) ~ sin ki cos K2 sin [ijjp — tl>i)\ An 

+ [cosKi sinK2 cos {ij}p — 1P2) — siuKi COSK2 cos (V'p — "01 )] Qn 
+ sin Ki sin K2 sin ('02 — "01 ) Ln , 
then 

I Si X S2 1 • f = cos Ki sin K2 sin (-^p + Xp — V'2) — sin ki cos K2 sin {-ipp + Xp ~ i^i) 
Si X §2 ) • Ln = sinKi sinK2 sin (■02 — "^i) • 
Thus we can rewrite Eq. ([57)1 in detail as 



^2^3 



3GL 



N 



■sin7 7 — ^ siuKi sinK2 sin (02 — "01) 

(1 — W2i'^ ) — — sinK2 cos (0 — 02) — (1 — wiv) — — siuKi cos (0 — ipi) 
X [cos Ki sinK2 sin (0 — 02) — siuKi COSK2 sin {ip — 0i)] , 



(37) 



(38) 



(39) 



(40) 



where tp — ipp + Xp- Again, it is manifest, that the relative angle of the spins stays constant for equal mass black 
holes, irrespective of their orientation. 

Starting from the above remark, Eq. (p6| and the estimates (JB35P we find 



O (7) = O (e) O (i/3-2») 0(r-i) 



(41) 



Thus the angle 7 changes faster than k^. 



Eq. (30) of ^ 



gives 0i = arctan 



C. Spin azimuthal angles tpi 



Si = sin Ki cos 0il + sin Ki sin -0^1X1+ cos KiL 



N 



(42) 



rh • Si ) / 1 1 • Si ) for the spin azimuthal angles, unless kj ~ 0, tt (the spins are aligned or 

antiahgned to the Newtonian orbital angular momentum) or ipi — 7r/2, 37r/2 (the projections of the spins in the plane 
of motion are perpendicular to the node line) . 

In the generic case the spin azimuthal angles evolve according to 



(1 + tan^ 0,) ■0, (i • Sij = (m - tan -0,1) • -n^i+ ( ^m - tani/>,— 1 j -Si 

As both 1 and m precesses about J7l, while Si about J7i, we find 

^Ai siuKi = ( cos0im — sin^J) • (ili — rJi,) x Si , 
or, by employing Eqs. (I^- (PT|) : 

0iSinKj= sin('0p -0i) An + cos(V'p -0i)QN • (Jli - Jl^) x Si 



(43) 



(44) 



(45) 



with the vector products Jli x Si and H^ x Si, given by Eqs. (jB39p and (jB6l) . respectively. We obtain 



■ipi = f Jli • Ln j + fjli • QnJ sin(V'p - ^i) - f Oi • AnJ cos(V'p - V'i) cot k^ 

LLT 

-as-^ [cot a sin {tpp + Xp) - cot Ki cos {xp + ipp - i^i)] ■ (46) 

With this we have completed the derivation of all required evolution equations. 
Starting from Eq. pS)) and the estimates (IB35[) we find 

©(tA,) =0(e)0(l,ry)0(T-i). (47) 

The change in the azimuthal angle of the spins is one PN order higher than the Keplerian orbital evolution. 



V. SPECIAL CONFIGURATIONS 

As a by-product of the calculations carried on in this paper we have recovered the known result that the plane 
of motion is changed only by perturbing forces pointing outside the plane of motion, thus by the SO, SS and QM 
perturbations. We have shown that the relative angle of the spins stays constant for equal mass black holes, irrespective 
of their orientation. We have also proven that unless the spins are perpendicular to the plane of motion (ni = 0), the 
polar spin angles will change under these perturbations. 

The nonprecessing (k^ — 0) and precessing (generic Ki) cases have been discussed sep arately in the literature (see 
Refs |47|1 and \4^ , respectively) in connection with the recoil of the final black hole [49| . From among the precessing 
cases the antialigned spin configuration with the spins laying in the orbital plane has received special attention, as 
numerical investigations have shown that it leads to the highest kick velocity. 

We have now the means to investigate such a configuration analytically. First we specialize to spins laying in the 
orbital plane, k^ = tt/2. After some algebra, Eq. ([M)) gives 



K. 



so 



j^SS ^ j^QM^^ 



K 



so 



G'^m?ri 
2c3r3 

_sin(V., + Xp-V^.) ^ ^4^,.._3 ^ 3) ^^ 



(48) 



1 






'Xp 



k=l 



2cos(-0p + Xp - V'fe) + 



A 



N 



Gnifi 



[2cos{i>p -Tpk) - 3 sin Xp sin (i/ip + Xp - V'fc)] 



K 



SS 



K 



QM 



1?' ^Xj [3 sin {2%l)p + 2xp - -0j - Vj) + sin (V-j - V'i)] , 
ZwiXi sin (2?/;p + 2xp - 2i})^) . 



All contributions Kf-*, Kf^, K^ are of the same order. In general the expression for ki does not vanish, not even 
in the special case of equal mass {v = 1), maximally spinning (xi — 1) black holes {wi = 1) on circular orbit {An = 0), 
when 



3-3— [2 sin (2V'p + 2xp - 2V'j) + 2 sin (27/'p + 2xp - ipi - i>j) + 3 sin (ipj - t/ji)] 



(49) 



Therefore in general a configuration with the spins in the plane of motion is not preserved. 

However in the special case Tpj — ipi + n and equal mass (v ~ 1), equal spin (x2 — Xi) black holes {wi ~ 1) we find 



as 

Hi An — 

J^i Qn = 

r2i ■ Ln ^- 



0, 



Xi cos {■)pp - ipi) 



- c3^3 Xism(^p-V^,) , 
7G 



T J cos a , 



(50) 
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such that according to Eq. (p4)) ki = 0.^ 

Then one has to check, whether the condition imposed on ijji is consistent with their evolution. With 03 = Eq. 
(|B6|) gives rii X Si = 0, while from Eq. (|B39|) we get 



7G 
Oi X Si = — — — -Jcosa 



such that Eq. (|46p simplifies to 



sin('(/'p - V'i) An + cos(?/'p - V'i) Qn , (51) 



7G 
^'= 2^2^:3^ cos a. (52) 

As the right-hand side does not depend on the index i, the imposed antialignment of the spins can be maintained 
over time. This is also evident from Eq. (pB)) . We have also checked that the constraints (46)-(47) of [4l[ are trivially 
obeyed. 

Therefore the special configuration of equal mass black holes with equal, but antialigned spins, both laying in the 
plane of motion is preserved by the conservative PN dynamics, with leading order SO, SS and QM contributions 
included. This stands as the main result of this section. 

Equation (48) of [41] allows us to rewrite 

7G 

i^i = r, 2 3 ^^ (^ + ^-P^ + *^2PAr) , (53) 

with the coefficients (given by Eqs. (39)- (40) of f4l|) specified for equal mass as 

1 /v\2 13Gm 



£PN 



^2PN 



3 /u\4 liGm fry 63Gm/w\2 / Gmy 



128 VcV 32 cV Vcj ^32"^Vc 

VI. CONCLUDING REMARKS 

In this paper we have established the conservative evolution equations of the two independent sets of variables 
characterizing a spinning compact binary during its inspiral, established in [4l|, with leading order SO, SS and QM 
contributions included. As the lengths J and Xi ^'^^ constants, this reduces to angular evolutions. The evolutions of 
the variables complementing the set (J, Xi)i the inclination a and the spin polar angles Hi were given as Eqs. (1^ and 
P4p . The evolution equations for the spin azimuthal angles "0^ (replacing Xi ^^ independent variables) were given by 
Eq. (1461) . These time derivatives (and all others computed throughout the paper) can be transformed to derivatives 
with respect to Xp by employing Eq. (J17p in the form 

^, = {—.-^A-U)^. (56) 

dt \fir-^ J dxp 

The true anomaly Xp becomes the only independent variable by employing the parametrization r (xp), Eqs. (|B20p - 

(mn). 

The system is closed by the evolution of the argument of the periastron i/jp given as Eq. (PTj) . the last two 
Eqs. ^ giving An and Ln, the analytical expression (jB17P - (jB19P of the perturbing acceleration components a^, 
the expressions ()B36p - (|B37p of the of the spin precessional angular velocity components I fti ■ An ) and ( Jlj • Qn ) , 

finally the vector products Jli x Si and fl^ x Si, given by Eqs. (|B39I) and (|B6p . respectively. 

Therefore we have derived a closed system of first order ordinary differential equations for the variables (a, Hi, tpi, 
ipp, An, Ln) evolving in terms of the true anomaly Xpi ready for numerical evolution. From this set (a, Ki, il)i) are 
independent variables characterizing the spinning binary configuration, while {ipp, An, Ln) characterize the orbit. 



^ The so contribution to ki vanishes, while the SS and QM contributions cancel. A glance at K~^ given by Eqs. II48I I shows that 
without imposing the black hole condition Wi = 1 the SS and QM contributions do not cancel, therefore the result does not hold for 
equal mass, identically spinning neutron stars. 
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In another way of counting, replacing {An, Lm) and their evolutions by the orbital elements (a^, e^) and Eqs. 
(P))- (IT(I1) . respectively; also including the evolution Eq. (P7)) for the longitude of the ascending node — </>„ we have 
obtained evolutions for (i) the orbital elements (a^, e,., a, ipp, —4>n) characterizing the perturbed Keplerian motion 
and for (ii) the spin angles (k^, ^i) characterizing the spin orientations with respect to this perturbed Keplerian orbit. 

As a by-product, we have proven that the relative angle of the spins stays constant for equal mass black holes, 
irrespective of their orientation. 

Also, unless the spins are perpendicular to the plane of motion, the polar spin angles change under the perturbations. 
There is one notable exception under this rule: the special configuration of equal mass black holes with equal, 
but antialigned spins, both laying in the plane of motion is preserved by the conservative dynamics. This is the 
configuration which led to maximal recoil found in numerical simulations [48|, and our investigations show that it is 
conserved during the inspiral to a 2PN accuracy, with leading order spin-orbit, spin-spin and mass quadrupole effects 
included. 
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Appendix A: Comparison of notations with related literature 



In this Appendix we compare the notations established in [4ll| and thoroughly employed in this paper with corre- 
sponding notations in the literature. 

First we establish the correspondence of the Euler angles (—(/)„, a, f/'p) employed in [41| and standard celestial 
mechanics angular orbital elements in Table HI The celestial mechanics angular orbital elements (J7,i,a;) are defined 
with respect to a reference plane and a reference direction contained within it, both inertial. The node line is defined 
as the intersection of the reference plane with the plane of motion; the angle span by it with the reference direction 
is the longitude of the ascending node Q; the relative angle of the two planes is the inclination l and the angle span 
by the ascending node with the direction of the periastron in the argument of the periastron uj. The Euler angles 
(— 0„, a, %l}p) employed in [4l| are defined similarly, but with respect to the inertial system /Ci with x and J standing 

as the X and z axes [any x ± J standing as the reference direction and the reference plane given by ( x, y = J x x ) 



TABLE I: Comparison of the notations in Paper I and standard celestial mechanics angular orbital elements. 
Rcf. [411 Celestial mechanics 



Euler angles ( — (^„ , q , f/ip ) 

True anomaly Xp 

Equation (14) of [41J <^„ = -q- 

Equation (15) of [41] "ijip + icv = 



{i'p^ 






+ ( 



angular orbital elements ($!,(,, a;) 
true anomaly v 



A . tan(Lij + u) 






The various systems of reference necessary for the description of the motion were also discussed in Refs. [43| and 
p3| . We establish the correspondence in Table HT] While in these papers a quasicircular orbit was assumed, the results 
of [4l| hold for generic orbits. A correspondence can be established as long as J can be viewed as an inertial axis. 

Finally we compare the notations of Ref. j4^] with the notations of [4l| in Table [ 



Appendix B: Decomposition of the acceleration and spin angular velocity vectors in the system /C4 during 

the inspiral 



In this Appendix we give the decomposition of the accelerations and of the precessional angular velocities of the 
spins in the system /C^. The ingredients we need are Eqs. (19)-(20) of 41] for the decomposition of the position and 
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TABLE II: Comparison of the notations in Refs. [43] and Q and [41 





Reference [43] (based 


on 


[42J) 


Reference [41] 


Orbit 


circular 








elliptical, An / 


Corresponding quantities 


{(jj, Muj) 








f Ln _ c^ ,3/2 c3,3/2^ 
[H;^ Gm^ ' G^ j 


Corresponding equation numbers 


(16), (18) 








(25) and (26), (15) 


Plane of motion 


orthonormal base (n, 


A) 




nonorthogonal base (7, f ) 


Orthonormal inertial source system 


(ef,e,^ef) 








7^, = (x,y,j),if J = ef 


Orthonormal basis in the plane of motion 


(ef,ef) 








("f, mV if J = ef 


Euler angles 


(*s,t,a) 








(V) = V'p + Xp, a, -1 - '/'n) , if J = ef 


Line of sight in the {x, z) plane 


© = arccos ( ef 


f,) 






e = arccos (j. at), y = £f 



TABLE III: Comparison of the notations in Refs. \4£ 

Reference [441 



and ^. 
Reference [41] 



Orbit 

Corresponding quantities 

Inertial axis 

Inertial orthonormal basis _L Jo (to J) 

Inertial orthonormal basis _L Ln (the Kl basis) 

Basis comoving with ^ 

Symmetric mass ratio 

Polar and azimuthal angles of Ln in Ki 

Phase 



circular 

{uJorb, V'^ = MiOorb) 

Jo 

(*, y) 

(xl, Vl) 

(n,A) 

u £ [0,0.25] 

(t,a) 

$(i) 



elliptical. An 7^ 



c-3/2 cf, 3/2' 



Gm ^ ' G ' 



J 

(*, y) 

f, Ljv X f 
r;G [0,0.25] 
(a, f - <^„) 



velocity vectors: 



f = cosxpAn + sinxj9QN , 

GmiL 
V = 

Ln 



-sinXpAN + cosxp 



A 



N 



GmfiJ 



Q 



N 



(Bl) 
(B2) 



In the system K-l the spin is given by Eq. (30) of [41']. By employing Eqs. (|20p - ([^T|) we rewrite it in the system ICa 
as 



We also need 



and 



Si = sin Ki 



cos {-tpp - ipt) An " sin {ipp - ipi) Q 



N 



cos K, L 



N 



f X Sk = cosKfc fsinXpAN -cosXpQnJ - sinKfcSin(V'p + Xp - V'fc)LN 



±^ Gmn 

V X Sk = — cos Kfe 

Ln 

Grrni 

— sm Kk 

Ln 



cosxp 



- I An +smXpQN 



GmiiJ 



A 



n 



cos (V'p + Xp - ^fc) + -7; cos {^p - tpk) 



-■N 



Jli X Si = a^-— 
Ln 



sm K, 



tana 



sin (f/'p - ipi) sin {ipp + Xp) + cos n^ sin Xp 



^N 



sm Kj 



cos (^/-p - i/'i) sin (fAp + Xp) - cos Ki cos Xp 



tana 
sin Ki sin (V'p + Xp - V'i) Ln| 



Q 



N 



(B3) 



(B4) 



(B5) 



(B6) 
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Eqs. (3)-(5) of [41^ give 



G 

c 



S, = —m^r/v^' ^Xi , 



G2 

C4 






(B7) 
(B8) 



1. Acceleration 



The general relativistic, SO, SS and QM contributions to the acceleration, with the SO part given in the Newton- 
Wigner-Pryce spin supplementary condition (NWP SSC) [23, Q, by employing Eqs. (|B7|) - (|B8I) are: 



with 



A 1 1 NWP 



- ass + sqa/ 



a.PN 



^2PN 



Gm 



Gm I r3 



Gm 3 

2(2 + 77) (1 + 3ry)w^ + -r^f^ 

r I 



f + 2(2 -ry)fv^ , 



- (12 + 297?) \^) + r; (3 - 4r;) 7;4 + ^,7 (1 - 377) r^ 

--77 (3 - 477) r'^v' - ^ (13 - 477) —X? - (2 + 257; + 2772) —7^2 
2 2 r r 



G771 



ri (15 + 47;) 7;2 - (4 + 4I77 + 87;^) 377 (3 + 277) r 



^NWP _ G m 1] 



^so 



S-SS 



^QM 



E (4-"-^ + 3) X. { 1^ (Ln • S,) f - (v X S,) + I (f X S,; 
2^^X1X2 { [(Si • S2) - 5 (f • Si) (f • S2)] f + (f • S2) Si + (f . Si) S2} , 



3G377l377 



3G^77l^77 



2c4r 



4^4 



fe=i '^ 



1-5 f-S 



f +2 f • Sk S 



k ^k 



(B9) 



(BIO) 



(Bll) 

(B12) 

(B13) 

(B14) 



After inserting Eqs. (jBl[) - (|B5p . the projections a^ = Aa • f (j) with f(i) = (An, Qn, Ln), they can be readily found. 
For explicit expressions we also need 



f • Sk = sin Kfc cos {tpp + Xp - V'fc) , 
81-82 = cos Ki cos K2 + sin Ki sin K2 cos (7/^2 — V'l) . 

The acceleration components are 



(B15) 
(B16) 



- ^PN ,,.2PN ,^,SO ,„SS ,OM 



ai — a^ + a 



a, + a^ 



a^ 



-,PN 



2PN 



Gr 



2(2 + T])-!^ - (1 + 3t])v^ + -f]r^ 
r 2 



, . . Gnifi , 

cosxp - 2{2 - T])r— — smxp 



Gm 
3 



3 / Gm \ 15 

- (12 + 2977) f — j + T] (3 - 4rj) v'^ + —T] (1 - 377) f-* 

77 Gm 2 ( 2\ ^"^ -2 



„77 (3 - 4r/) T^^w^ _ ^ (13 _ 4^) „2 _ (2 + 25ry + 277^) r 

2 2 r ^ ^ 



Gtti 



77 (15 + 4r/) i;2 - (4 + 41r/ + 877^) 37? (3 + 2r/) ?^' 



r 
Gm[if 



cosxp 
sinxp > 



CO _ G^7n^77 



\ 2^r ^ T 



N 






3r . 
'Xp + TTSmXp 



2Lm 
2 



k=l 



af' = 



3(^3_3 



mi] 



4 4 XiX2i cosKi COSK2COSXP + sinKi sinK2 

X [[cos {lp2 - -01) - 5 COS (f/ip + Xp - 1pl)cOS {-Ipp +Xp - "02)] COS Xp 

+ COS (7/;p + Xp - V'2) cos (f/ip - f/ii) + cos {ipp + Xp - "^i) cos {ipp - 1P2)] | , 



,QM 



3G3 



777, 77 



2c4r4 



^ 7«fe7^^''"^Xfe|cosXp - sin^ Kfc cos (fAp + Xp - V"*:) 



fe=i 



X [5 cos Xp cos {ipp + Xp - V'fc) -2 cos (7/;p - tpk)]j , 
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(B17) 



« — «-P^ ^ «2PAr , so , „SS I „QM 
tt2 — ^2 I 0^2 \ ^2 ' ^2 ' ^' 



PW 



2PN 
0-2 



Gm 



2 ' 

2(2 + 7;)— -(l + 37;)7;2 + ^777^2 
r 2 



smxp + 2(2 - ijjr— cosxp ' 



^N 



Gm 



^ (12 + 29r/) ( ^ j + 77 (3 - 477) v'^ + ^77 (1 - St]) r'^ 



-r) (3 - 477) r^v^ - ^ (13 - 477) -^v^ - (2 + 2577 + 2ri^) -^r^ 
2 2 r r 



smXp 



Gm 



ri (15 + 477) w^ _ (4 + 4i,y + 877^) ^:^^ - 377 (3 + 277) f 



Gmiir 



,so 



ss 
"2 



G^m.'^T] 



3G m, rj 



3Ljv Gm,^\ 3r 



2Ln 
"^ (41^^''-^ + 3) Xk cos Kk 



Xp 



fe=i 



'-{' 



4 ^^ — X1X2 ] cos Ki COS K2 sm Xp + sm ki siii K2 



X [[cos (V'2 - t/^i) - 5cos (7/;p + Xp - V^i) cos (t/^p + Xp " V'2)] sinxp 

- COS {tpp + Xp - "02) sin (7/;p - ipi) - cos (7/;p + Xp - "^i) sin (V^p - 7/^2)] | , 



,QM 



3G%i^ 
2c47 



^M;fcp2''"^Xfe|sinXp - sin^ Kfe cos (V-p + Xp - "V/c) 



fc=i 



[5 sinXp cos (t/^p + Xp - "Vfe) +2 sin {tpp - ipk)]} , 



Gmn J 



A 



N 



GmiiJ 



(B18) 
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and 



n - n^O , ^SS . OM 
0-3 — 0^3 + ^3 + "^3 I 

G2 2 '^ 
u^ _ m rj 'sr^ fA,.2k-3 



(B19) 



^(4i/2'=~3 + 3)xfcSinKfe 



■,SS 



OM _ 



fe=l 
3G3^3 



N 



A 



N 



COS {ipp + Xp- ^^k) + T^ COS (-0P - f/ife) 



if 



sin (V'p + Xp - fAfc) ^ , 



4 4 -XiX2 [cos Ki sin K2 cos (-0P + Xp - '02) + cos K2 sin ki cos (V'p + Xp - ■0i)] 

oG ?7l Tj ^— \ 2fe — 3 2 ■ n / ; i \ 

' 2 4r4 Z^^fc'^ XfcSin2KfeC0s(V'p + Xp- V"*;) ■ 
"^ ^ fc=i 



In the above expressions we still need to employ Eqs. (20)-(21) and (23) of [41 1 



r2 



/i [Gmfji + ^jv cos Xp) 
^^smxp , 
2 {Gra^) + A^ + 2Gm^AN cos Xp 



r = -^smxp 



(B20) 
(B21) 
(B22) 



in order to rewrite r, f and u^ in terms of the chosen dynamical variables. 

Also, as Ln is not among the chosen independent variables, we need to express it in terms of them. For this, first 
we give the SO part of the orbital angular momentum in the NWP SSC: 



Gn 



T NWP 

^«o - 2? 



^ (4 + 3^.3-2*) S, 



f X ( f X Si 



fc=i 



G 



2,yy^6 



Ac^j 



if ^ (4j/2*"^ + 3) Xi (sin m [cos (2xp + i^p - ^i) - cos (i/'p - V'^)] ■ 



^N 



+ sinKi [sin(2xp + V'p - i'i) + sin(V'p — tpi)] Qn - 2cosKiLN \ 



As expected its order is 



'T NWP 

Ln 



= 0[e'"^)0{i^)0{l,u- ")x^ 



2i-3\ 



(B23) 



(B24) 



The total angular momentum J = L + S1+S2 in the system Ka then becomes, using Eqs. (|B3p and (|B23p : 

,2 2 



J J = ?7 2_^ Xi sm Ki 



Gm 



v^' 3 cos {iPp - iPi) - -^-^T] {iiy^' 3 + 3) [cos i^pp - V'j) - cos (2xp + V-p - A)] 
v^'-^ sin (V'p - V^») - ^ry {Av^''^ + 3) [sin (V^p - i^i) + sin (2xp + ipp - V'i)] 



Gm2 



2 r 



+ ILn{1 + epN + e2P7v) H rj ^ 



Gm 



-2*-3 ^"'.,,(4,y2»-3+3) 



2c2? 



An 
Qn 

Xi cos Ki S Ln , 



(B25) 



Here epjv and e2PN are given by Eqs. (39)-(40) of [4l|. The projections along the basis vectors 1, rh, Ln of the /Cl 



system are 



= ^ Xi sin K, 



Gm 



i ipi + -^-?7 (4j/2' ^ + 3) [cos (2xp + 2-0P - Vi) - cos tp^ 



Ac^ 



cJ sin a 
Gm?"q 



E 



= > Xi sm K., 



G 



m 



v^^'^ sin V'i + —^V (41^^' ^ + 3) [sin (2xp + 2?/^^ - i/-*) - sini/'* 



G 



4c^r 
2 2 



T /1 N <-^"* V^ 

J COS a = Ljv (1 + tpN + e2PAr j H r; _> 



4 = 1 



Gm 
2^ 



7] Uv 



2J-3 



Xi COS Ki 



The last equation enables us to express Lj^ to 2PN accuracy in terms of the chosen independent variables: 

Ln = J {^- <^PN - (■2PN + CpAr) COS a 



Grn^ 



i=l 



(1 - eppf) v 



,2i-3 



Gm 
2^ 



ri (4i/2'-3 + 3) 



Xi COS Ki 



We also give the series expansion of its reciprocal: 
1 _ 1 + epTv + e2P7V , f Gm^ 



+ 



J cos a 
/ Gm^ 



V cJ cos a I J cos a 



V 



Gm 

(1 + epN) Xv - T^n \^Xv + 3x+) 



?7" 



Gm, 



(1 + epjv) Xz^ t"*^ (^^'^ + 3^+) 



\cJ cos a/ J cos a 

f Gm^ \' 773 3 ^ Gm2 ^" -^ 



Xi^ 



-X. + 



?? 4 



\cJ cos a/ J cos a "^ \cJ cos a J J cos a 

where we employed the notations 

2 

X+ = ^Xi cos Ki = Xl cos Ki + X2 cos K2 



'i=l 

2 



Xl- 



/^ '^^^ ^XiCOSKi = J^ "'^Xl COSKi + i^X2 COSK2 . 



j=l 
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(B26) 
(B27) 
(B28) 



(B29) 



(B30) 



(B31) 



Note that the 2PN contribution of 1/Ln is rather messy (fourth rank in the spins), nevertheless for our purposes we 
need it only to IPN accuracy (it enters only in FN terms or higher, and the desired accuracy is 2PN). 

We also give here the detailed expression in terms of orbital elements of epjv, which is necessary at this accuracy: 

1-377 /w\2 Gm 

epN = — 7, — (-) +3 + 77)^- 
2 Vc/ c'r 

^2 



(7 - 7/) {Gm^i) + (1 - 377) A% + A{2-ri) Gm^xAjq cosxp 



1c Lj^ 



{Gm[i) 



[(1 - 37/) e^ + 4 (2 - 77) er cos Xp + (7 - ??)] 



(B32) 



2. Spin angular velocity 

The spin undergoes a pure precession, therefore its magnitude is unchanged, while its direction changes as 

d 



dt 



Si = fliX Si 



(B33) 



These are the equations in the NWP SSC corresponding to Eqs. (46)-(48) of 41], which were written in the covariant SSC. 
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where, after employing Eqs. (jB7|) - (IB8|) . (jBip . (|B3[) . and (JB15I) in Eqs. (56) of [4i| the angular velocity vector is found 
as 



fti 



n 



so 



n 



ss 



fl 



QM 



^ 0(4 + 3^^ ^ 



2c2^3 



fl 



SS 



2c3r3 ^^ 



TV 

sin^ 



N 



Kj I [3 COS {ipp + Xp - ip] ) cos Xp- cos (?Ap - i/'j )] An 
+ [Scos (V'p + Xp- V'j) sinxp+sin (-(Ap - -0^)] Qn| - cosk^Ln , 



m 



QM 



Gm^rj 



iwtXt sinKi cos (-0P + Xp - V'l) f cosxpAn + sinxpQN 



2c3r3 
with j ^ i ■ Their PN order is 

0(nf°) = C'(e)e»(l,i^3-2')C(T"^) , 
0(f2P) = O (£3/2) 0(77)0 (1^2^-3) X,C3(T-i) 

0{n9''^ = o(e'''')o{^)w,x^O{T-^). 

The projections employed in the main text are 

Oi-AN = 33 {v^^^^Xj sin Kj [3 cos (-0p + 2xp - V'j) + cos (V'p - V'j)] 
+3wiXi sin Ki [cos (V'p + 2xp - i>i) + cos (V^p - Vi)]} 7 

2^3^3 {^'^■'"^Xi sin Kj [3 sin (V'p + 2xp - "^j) - sin (V'p - Vj)] 
+3«7iXj sin Ki [sin (V'p + 2xp - ^i) - sin (V^p - -0^)]} , 



rji Qn = 



rJi-LN = — „ , ^ ^ J cos a- „ ^ ^ [(4t/2» 3 _,_ 3) ^^ (^og ^^ _,_ j^2j 3 (^5 ,,.3^3 )Xj' 



(B34) 



(B35) 



We also need 



2c2r3 



Oi X Si = 



2c3j,3 



(B36) 

(B37) 
(B38) 



rJi ■ Ln ) sin Ki sin (f/^p — 0^) + I Oj • Qn I cos k., 
rJi • Ln ) sin Hi cos (0p — ipi) — [fli- An ) cos k. 

Uli ■ An j sin (0p - t/;,;) + f Jli • QnJ cos {ipp - ipi) 



An 
Qn 



sm K. 



-'N 



(B39) 
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